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Abstract 

We calculate the leading quantum and semi-classical corrections to the New- 
tonian potential energy of two widely separated static masses. In this large- 
distance, static limit, the quantum behaviour of the sources does not con- 
tribute to the quantum corrections of the potential. These arise exclusively 
from the propagation of massless degrees of freedom. Our one-loop result is 
based on Modanese's formulation and is in disagreement with Donoghue's re- 
cent calculation. Also, we compare and contrast the structural similarities of 
our approach to scattering at ultra- high energy and large impact parameter. 
We connect our approach to results from string perturbation theory. 
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I. INTRODUCTION 



Quantum gravity is a subject beset with calculational difficulties, most predominantly 
due to ultraviolet divergences that are incumbent in any 4-D theory with a coupling constant 
with dimensions of length. These difficulties are well documented in the literature and have 
not allowed a consistent formulation of the theory at short distances . The hope is that 
superstring theories will provide us with a robust definition of the theory in that regime. 
However, as Donoghue [Q] recently pointed out, one can analyze quantized general relativity 
within the computational framework of effective field theory. He argued and showed that the 
unknown ultraviolet structure of the theory is irrelevant for the purposes of obtaining well- 
defined leading-order in h quantum corrections to some quantities involving gravitational 
effects at large distances. It is possible to formulate this program in a general way armed with 
techniques based upon general analyticity properties and gravitational Ward-Slavnov- Taylor 
identities. As we will discuss, this strategy is also implicit in the analysis of the scattering 
problem at ultra-high energy and large impact parameter considered by a number of authors 
in the past few years. The problem of long range forces is primarily a consequence 
of soft virtual gravitons and inherently of an infrared nature. The infrared region with 
respect to radiation of soft gravitons and the Block-Nordsieck [^] cancellation of infrared 
singularities was communicated long ago by S. Weinberg ^ and B. DeWitt 0. 

Donoghue employed this method to obtain the leading corrections to the Newtonian 
potential. Our starting point is similar in spirit but quite different in implementation. We 
present and apply a general formalism presented recently by Modanese |^ which we believe 
is better suited for the study of the simplest quantity in the long distance regime, namely 
the static potential. We obtain a different result than Donoghue, and we discuss the reasons 
for the discrepancy in detail. In particular, we compute the quantum corrections in the 
gravitational or more generally the zero-mass sector of the theory in a manner that can be 
applied to either the static or ultra-relativistic problems. This is in parallel with studies 



of non-Abelian gauge theories. [|T0]-[T3|. The approach of Modanese is presented in terms 



2 



of the Euclidean functional integral and parallels the Wilson loop of gauge theory. The 
formalism is quite general and allows a starting point for perturbative investigations of the 
long distance regime of this work as well as non-perturbative numerical simulations. On the 
other hand, the application to the relativistic scattering domain has been communicated in 
the literature under the names of relativistic eikonal, and semi-classical (WKB) techniques by 
many authors fl^-Q . The two problems are very closely related and a general result can be 



derived by functional integral techniques which can be characterized by the exponentiation 
of certain connected Green's functions. The connected nature of the exponentiation of the 
Green's functions separates the potential from its iterations in a clear manner. 

The actual work in such calculations is two fold. Firstly, one identifies a formulation that 
allows a convenient isolation of the effective theory in the relevant kinematic domains. In the 
static limit the relevant effective theory is a 3-D, dimensionally reduced quantum gravity, 
while in the ultrarelativistic domain, the effective theory is 2-D. Secondly, the isolation 
of contributions that survive as inverse powers at large distances, i.e. as R~"' for n > 1 
correspond to singular contributions at small momentum transfer i.e. analytic in the cut 

plane with branch cut singularities at zero momentum transfer. The Fourier analysis of 
the Newtonian potential (n = 1) is already singular at small momentum transfer with the 
well known pole in Q^. The specific power and nature of the singularities depends upon the 
number of graviton interactions and other details that will be elucidated in the following 
sections. 

Contributions analytic at small momentum transfers (in the sense of old S-matrix theory) 
correspond to forces of finite range of the form exp— due to the Riemann-Lebesgue 
Lemma, and permit a Taylor series expansion at small momentum transfer. It is tempting 
to say that modern superstring theory will relate the range or the effective length parameter 
to the string tension in a way that all such contributions can be rendered finite. However, 
tempting such a correspondence is, its proof certainly lies outside the scope of this investi- 
gation. Nevertheless, we will expound on this idea using known one-loop results from string 
theory in a later section of this communication. The field theory limit with zero mass par- 



tides in loops arises from string interactions (higher genus Riemann surfaces), as required 
by unitarity. This hmit is achieved at zero string tension, with polynomial corrections with 
finite coefficients which characterize the surviving string interactions at low energy. It is 
only the effective field theoretic (or S-matrix) contributions with zero mass intermediate 
states that reduce to non-polynomial interactions at low energies. 

The organization of this communication is as follows: In Section II, the underlying 
technique pioneered by Feynman 0] and developed by DeWitt [§] and Mandelstam [l^ is 



briefly reviewed. In this approach the metric is treated as a tensor fleld on a flat space-time 
background in parallel with any other quantum fleld. Path integral quantization is used as 
a convenient framework to introduce external sources to formulate the static potential (or 
scattering amplitude). The recent work of Modanese is presented, as well as mention of 
the older work of many authors on the high energy scattering. The relevant dimensionally 
reduced effective fleld theories are identifled. A discussion of the results in terms of connected 
graphs of the gravitational and more generally the zero mass sector is presented. Section II 
also deals with the calculational aspects of this work, namely the relevance of the different 
distance scales and the isolation of the power behaved contributions at large distances. 
Calculations in the ultra-relativistic regime, relevant to the scattering of constituents at 
Planckian energies are reviewed but not presented in any great detail, since these have been 
dealt with in the past few years by a CERN group ||19|| . 



Following the explicit calculations, a general anatomy of a proof is presented in Section 
III to indicate how one would make a general proof based upon the plausible scenario that 
the small distance effects do not survive at large distances. Here recent calculations of higher 
order (higher genus surfaces) contributions in string theory are quite necessary. 

Finally, we consider the pitfalls and limitations of such an investigation and compare and 
contrast gravity with other gauge theories, where long range forces play out in different ways. 
In abelian gauge theories the flnite fermion mass screens all contributions other than the 
Coulomb potential between charged particles in the deconflned phase. There are also non- 
trivial, inverse-power behaved van der Waals forces between neutral objects resulting from 



effective non-renormalizable interactions at low energies . While the techniques in gravi- 
tational theory are similar to the analysis of molecular forces, the physics is vastly different. 
In non-abelian theories, confinement and spontaneous symmetry breakdown are plausible 
dynamical mechanisms that render long range inverse power behaved forces inoperable. An 
appendix deals with some technical questions on the details of the calculations. 



II. FORMULATION 
A. Static Limit 

The general expression for the static potential energy for the gravitational interactions of 
two static massive particles requires some care and has been given recently by G. Modanese 
0. His work builds upon some expositions of other authors in both the static limit and 
the ultra-relativistic limits. Standard conventions will be followed with a fiat space metric 
signature given by rj^i, = diag(l, — 1, — 1, — 1): we will work in M4 (4-D Minkowski space) 
instead of Euclidean space, since we will be interested in the perturbative definition of the 
theory with no applications to numerical simulations, or issues of instability. The Einstein- 
Hilbert action is given by 

-^EH = - / dx ^)R{x) , (2.1) 

where the integral is over the full four dimensional space time, g{x) = —det{gfj_^) and R 
is the curvature scalar. It suffices to consider the action of two constituents, since we are 
considering the two body potential. The coupling of the sources to the metric is via 

Sm = -Ml J dsi -M2J ds2 , (2.2) 

where ds = {g^y{x{T))x^x'^)'^dT is the line element along the trajectory of each object of 
mass Ml and M2 respectively, r is the proper time measured along the trajectory, and x 
is the covariant four-velocity. The generating functional for the gravitational interactions 
between two such massive particles is given by the functional integral: 
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\\Dg\ exp^^EH 

The integration is over the metric g with some suitable gauge fixing for the perturbative 
definition of the Feynman diagrams of the theory. We seek the potential energy of two static 
sources situated at 

x;' = (ti,0) x^ = (t2,^), -T/2<U<T/2 (2.4) 

These trajectories correspond to two static particles with coordinate separation R and four- 
velocities = 6q, with no relative motion. The limitations of the restriction to static 
sources within the context of the full dynamics is discussed below. 

In the perturbative definition of the functional integral Eq. (2.3), the metric is expanded 
around the fiat space metric r/^^^ in the usual fashion of the Feynman-DeWitt method of 
quantization of general relativity, namely 

gf,,y{x) = rif,y + Kh^^{x) , (2.5) 

where = 32ttG. Following Symanzik [^, we shall identify the lowest energy of the system 



(the gravitational potential energy, in the static case at hand), with the limit 

V{R)=\\m ?^lnZM, (2.6) 

in view of the identification of Zu with the vacuum-to- vacuum amplitude in the presence of 
external sources, given by Sm, adiabatically switched on at —T/2 and off at T/2, 

Zm = (0|e-*^^/^|0) = ^ e-*^"^/^|(0|n) p , (2.7) 

n 

for a very large time T. It is clear that Eq. (2.6) above is meaningless as it stands. We give 
it meaning by endowing the large Minkowski time interval T with a small negative phase, or 
by keeping track of the ie's in the perturbative propagators. This particular definition of the 
ground state energy was first advocated by Symanzik and proven in perturbation theory for 
the case of a source coupled linearly to the quantum field. It has been used in many other 



field theory contexts. In the apphcation to the static potential problem, this formulation 
clearly differentiates the static potential from its iterations, the higher order terms in T of 
Eqs. (2.6) and (2.7). We will assume that this definition holds even in the case presence of 
sources coupled to higher-order local monomials of the field. 

If g is expanded as in Eq. (2.5) in Eq. (2.3), we obtain (to lowest order in n) the customary 
result, 

Zm = (exp J dx h^,{x)T^''{x)), (2.8) 

where T^jy is the classical energy-momentum tensor of the source, and averages with respect 
to the normalized gravitational functional integral are defined in the usual way, 

J [Dh\ exp ^Seh 

A trivial time dependence proportional to exp—i{Mi + M2)T/h may be absorbed into the 
normalization of the functional integral. 

In the case of a pair of static heavy particles of masses Mi and M2 with trajectories 
given by Eq. (2.5), we have 

= -Ml j dti^l + Khoo{x{ti)) -M2I dt2^1 + Khoixit2)), (2.10) 

and 

T'^^x) = M,6^6',6^^\x - x{n)) + M26^6'q6^^\x - xir^)) . (2.11) 

We close this subsection with some remarks. In general the definition of the potential 
energy Eq. (2.6) and Eq. (2.9) appear to be frame-dependent and path-dependent in any 
field theory including gauge theory. The analogue in gauge theory, the Wilson loop, is 
gauge invariant but path-dependent unless the field strength tensor vanishes identically^. 
While the analogous expression in general relativity does not appear generally covariant 



^It is possible to construct a reparametrization-invariant "Wilson- loop like" quantity (see 
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(gauge invariant), the results Eqs. (2.3) and (2.10) are still a fully covariant definition of 
the interaction in the static limit much in the same way that the Coulomb interaction is 
a necessary ingredient in the choice of Coulomb gauge in quantum electrodynamics for a 
consistent gauge invariant dynamics. Of course it will be necessary to consider kinetic energy 
and the full dynamics to establish the general covariance in a completely convincing manner. 

The particular choice of trajectories Eq. (2.4) insures that the static potential alone is 
being calculated, and that no terms involving the kinetic energy or relativistic corrections 
are being included. At the level of {v/cY and V ■ {v/cY relativistic corrections are bound 
to become important in the determination of motion in accord with the post Newtonian 
approximation of classical general relativity, where both terms of order GM/Rc^ and {v/cf 
are retained. The virial theorem of classical mechanics relates kinetic energy Mv^ to poten- 
tial energy V. Therefore, kinetic energy will have to be included to order {v/cY in order 
to have a consistent dynamics for the purpose of computing the subsequent motion of this 
pair of gravitating objects. In our work Eq. (2.6) and the incumbent static sourses Eq. (2.4) 
form a convenient starting point for the computation of the static potential. The last remark 
underscores the statement above concerning gauge invariance (general covariance). 

Notice that in the weak couphng hmit Eq. (2.8), the lowest-order energy is invariant 
under infinitesimal coordinate transformations i.e. 



for a conserved energy momentum tensor in Eq. (2.8). The complete static problem also has 
a partial covariance with respect to time independent coordinate transformations that are 
in the little group of gQo{xi{t)), where Xi{t) are the trajectories, Eq. (2.4). 

It is also worth mentioning that neither Eqs. (2.6) or (2.8) is a loop integral, since the 
notion of oppositely charged sources does not exist for a tensor theory of gravity. This 
property is essential in obtaining the Wilson loop for the static energy in gauge theories. 
One may also consider the possible physical meaning of a gravitational loop functional or 
holonomy. 
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U{C) = -4 + TrP(exp i <f T^dx") , (2.13) 
where one views the connection coefficients matrix elements of a matrix F^. This 



issue was investigated by Modanese . He showed that the expectation value of the loop, 
calculated perturbatively around a flat background, does not give rise to long-range interac- 
tions, and so cannot be interpreted as representing a gravitational potential. The physical 
interpretation of such a result can be understood in terms of the functional integral. To 
leading order in h, the vanishing of the holonomy implies that the weak field configurations 
which contribute to the functional integral have zero curvature (to that order). An intuitive 
way to understand this discrepancy between gravity and gauge theories is by noting that 
the coupling of the vector field to matter is in terms of / df'xA^j^. For a pair of oppositely 
charged static sources, this coupling reduces to the Wilson- loop exponent. For gravity, on 
the other hand, the coupling to a source is non-linear in the metric (see Eq. (2.10)). Even 
the perturbative expansion involves a coupling of h^^, with the energy-momentum tensor T^y 
(see Eq. (2.8)). The affine connection coefficients simply do not arise in such an expansion. 

We proceed next to the calculations of the ffist-order corrections in G, both semi-classical 
and quantum mechanical, to the classical Newtonian potential energy, using the formalism 
developed in the previous sections. The choice of trajectories in Eq. 2.4 ensures that the 
static limit is retained as opposed to calculating any effects associated with relativistic 
corrections. In the following the harmonic or de Bonder gauge is maintained, 

5|.(/i^'^-^^^'^A) = 0. (2.14) 

In the harmonic gauge the graviton propagator is simply 

D^upa{x) = P^upa^ix) = ^{il^prj^^ + rj^^rj^p - ri^yr]p^)A{x) , (2.15) 



where 
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is the free massless propagator. The choice of gauge restricts the form of Green's functions 
and the perturbative Feynman rules as expected. An appendix deals with these and related 
issues in some detail. It is also assumed that the cosmological constant is zero. This is very 
important, a non-zero cosmological constant is an obstruction to the long distance behavior 
of Green's functions and their incumbent Ward-Slavnov- Taylor identities. 

Next we apply the strategy given in Eqs. (2.3) and (2.10) to calculate the potential to 
order G^. The basic strategy is to expand the matter action Eq. (2.10) to order G (n^), 

Sm = -(Ml + M,)T + J2 ( J dUhooixiU)) + dUh,,{x{U)f + . . . j (2.17) 

2 

= _(Mi + M2)T + T.{'J^J dtihoo{x{ti)) + KiJ dtihoo{x{ti)f + . . }j (2.18) 

and expand the partition function in powers of the sources Jj and Ki. Next connected 
graphs in the gravitational sector are retained. The graphs are connected to sources which 
are vertices for graviton lines. The complete set of graphs contributing to the static potential 
to order k'^ (G^) is shown in Figs. 1 to 3. The time coordinates of the sources are to be 
integrated from — T/2 to T/2 and we are left with an effective 3-D theory. The sources of 
mass Ml and M2 are coupled to the gravitational sector at positions and R. The potential 
or static energy can be read off from Eq. (2.6) in the large T limit. 

The connected nature of the Green's functions of hoQ{x) and local monomials thereof is 
extremely important and differentiates the potential from its iterations, which are terms of 
higher order in T in the expansion of the exponential in Eqs. (2.17), (2.18). 

The contribution to lowest order in G can be computed from the single graviton exchange 



graph in Fig. 1 [^. The result is straightforward and familiar 



Expressing the integrand as a contour integral and integrating gives (in the limit of large T) 

/i = ^T— (2.20) 

which, in view of Eq. (2.6), yields the Newtonian potential energy 
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V(R)^-'fM^. (,21) 

The corrections to order are of two classes: classical general relativistic (independent of 
h) and quantum mechanical of order h. They are enumerated in Fig. 2 for the semi-classical 
contributions and Fig. 3 for the quantum contributions. It is easy to understand that a 
perturbative expansion of the functional integral in gravity gives us classical contributions, 
both formally and intuitively. The formal reason is that the expansion includes tree graphs 
connected to an arbitrarily high number of external classical sources. Boulware and Deser 
2^ showed that the tree graphs reproduce classical general relativity (see also [p5| ,p6|] ) . 



The intuitive reason is that in gravity there are two dimensionless parameters that lend 
themselves to a perturbative expansion, one which involves % and is the ratio of the Planck 
length with the separation of constituents, and one which is independent of h and is the 
ratio of the Schwarzchild radius of a mass with the separation of constituents. 

The quantum contributions involve closed loops of gravitons and are of order h. Fig. 3a 
is the vacuum polarization and receives contribution from the Fadeev-Popov-FejTiman ghost 
which is a vector particle in gravity; this contribution is well known in the literature and 
has been calculated long ago. 

The calculations are discussed in detail in an appendix. We note the results here. The 
corrections to the Newtonian potential are denoted by 5V and the results are as follows. 



GM1M2 ( G{M^ + M2) 1, 
R \ Rc^ ^ ^ 2' 
GM1M2 ( G{Mi + M2) 



(2.23) 



R V 2i?c2 

from the graphs of Figs. 2a and 2b respectively. 
The quantum correction is of the form 



GM1M2 ( Gh A3 1 5A , , 

GMM 17 an ^^^^j 



R 20n R^(? 
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The detailed numbers come from Figs. 3a, 3b, and 3c respectively. In all of these contribu- 
tions, there is no difficulty in extracting the large time behavior in the leading static limit 
for large masses neglecting recoil in Eq. (2.6). Putting everything together gives to order 



All other contributions from skeleton graphs to this order are of short range and have 
support in a region where we do not believe the calculation in any case. 

We close this section with some remarks about the regime of validity of our results. 

• There is a hierarchy of distance scales in the problem, namely the Compton wavelength 
of the sources h/MiC, the Planck length {Gh/c^y^'^, and the Schwarzchild radii of the 
sources, 2GMi/(?. These scales are clearly not independent, since the Planck length is 
the geometric mean of the other two length scales. The corrections to the potential are 
to be trusted for distances large compared to all of these distance scales. This follows 
from the explicit expression of the corrections, but most importantly it is implicit 
in our restriction of staticity. Let us explain this more fully with a familiar Feyman 
diagramatic approach. Consider a particle interacting with a massive center through 
the exchange of a scalar particle (for convenience, not necessity) . If the momentum of 
the exchanged scalar is q and the initial and final momenta of the scattered particle 
are p and p', then the scattering amplitude is proportional to the sum of the direct 
and crossed propagators, namely 

^ + ^ . (2.27) 

(p + qY — m? + ie (p' — q^ — m? + ie 

In the limit of large mass and assuming on-shell scattering, po ~ Po ~ and the 
above sum becomes 




(2.26) 



2e 



12 oc 5(go) ■ 



(2.28) 



2mgo + —2m,qo + ie e 



'0 
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Numerator factors of go have been neglected in addition to the customary hnearization 
of propagators (Recall that spin 1 and spin two exchang will involve polynomials in 
p + p'. In this limit, one can sum a// the possible exchange contributions to scattering, 
in view of the following identity ||2^ 



X 



■ ■ ■ X ■ ■ — = (27r)"5(gio) x ■ ■ ■ x 5(g„o) , 

gio H \- In-ifi + «e 

where the delta function on the left hand-side is the overall energy conserving 6- 
function SiJ^Qio + Po — Po) ~ ^(J^Qio)- This identity shows that the sum over the 
dynamical exchanges is equal to a single graph in which the scattered particle field 
is replaced by a static pointlike classical source j{x) = 6{x) = F.T.[27iS{qo)]. This 
argument shows that we can neglect the dynamical, quantum behavior of the field 
which undergoes the scattering, provided that we are at a distance larger than its 
Compton wavelength and larger than distance scales where non-perturbative effects 
become important. The contribution of powers of go in the numerator will enter at 
order or higher in accordance with expectations on the low energy behavior of 
quantum electrodynamics and quantum gravity. It is also clear that these tecnniques 
can be generalized to establish a result that allows the factorization and the incumbent 
exponentiation for any finite correlation of exchanged quanta coupled to the source. 
This argument, which is not presented in detail here, is useful in establishing the 
general result in terms of connected Green's functions Eqs. (2.10, 2.11, 2.17, 2.18) for 
higher polynomials in /loo- 

The above argument implies that the corrections to the potential are valid for distances 
larger than the Compton wavelength of the source, if the mass of the source is less 
than the Planck mass {10^^ g) (with recoil corrections of order /Mi). This includes 
all of the elementary particles of the standard model including the hypothetical vector 
bosons that mediate the baryon number non-conserving interactions of Grand Unified 
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Theories and other massive particles of broken supersymmetric and technicolor theo- 
ries. For sources with macroscopic masses, on the other hand, the distance should be 
greater than the Schwarschild radius of the source. 

While the main results are similar in spirit to the work of Donoghue 0, the details 
and results differ. In our work the coupling to static sources Eqs. (2.6) and (2.10) 
is maintained throughout the entire calculation. In the calculation of Donoghue, the 
scattering amplitude with recoil and relativistic corrections in the coupling to matter 
are retained. The potential was defined as the non-relativistic limit of the one particle 
reducible graphs of the t-channel, a definition substantially different than the one 
presented here. An important check on the classical correction to the Newtonian 
potential obtained here is the complete agreement with the post-Newtonian potential 
energy of a test particle of mass M2 in an external field produced by Mi, in the limit 
Ml ^ M2 and where M2 has zero velocity, namely V = ^ + 1/2$^, where $ is the 
Newtonian potential . Our result for the potential energy is also consistent with the 
expansion of the 00 component of the metric in a Schwarzschild geometry in isotropic 
coordinates, 

2GM ( GM\ 

with the proper post-Newtonian identification ^fQo = 1 + 2$ + 2$^. 

The quantum corrections in Eq. (2.26) have a definite sign. The physical origin of this 
sign is due to positivity of the Hilbert space metric over physical states and is discussed 
in more detail in the appendix. Notice that the Newtonian limit is approached from 
above and the attractive interaction increases as the distance decreases. This is in 
agreement with the analogous situation in quantum electrodynamics with the form of 
the Uehling potential ||2^ . 



A final word of caution is in order. It is only at the one-loop level that the finite long 
distance behavior can be separated from the divergent small distance counterterms 
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in an unambiguous way, at least with the field theory techniques presented so far. 
Non-leading long distance quantum corrections at two loops and higher will receive 
contributions from divergent higher dimension polynomial counterterms. In order 
to completely understand the quantum contributions to long range forces in higher 
orders, one would need an improved definition of the theory. String theory is perhaps 
a suitable candidate. We return to this question in section III. 

B. Ultra- Relativistic Limit 

The study of the elastic forward scattering of two scalar particles at high energies greater 
than the Planck mass, s ^ M|j ^ t, with s = {pi + ^2)^ and t = {pi — ^2)^, has received 
attention in recent years because of its relation to superstring theory. 

One of the authors and M. Soldate investigated the problem ||^ because of the unitarity 
difficulty due to the energy growth of the scattering amplitude in any theory of gravity in 
near the forward scattering, 't Hooft and others were interested in the problem because 
of its possible relation to black hole dynamics, to string theory, and to other field theory 
contexts. The scattering can be treated nicely by semi-classical methods which go under the 
name of eikonal approximation or WKB and are well developed in the literature. The basic 
strategy is to retain the leading terms in the coupling of the external particles to gravitons 
at high energies and linearize the propagators. Probably the most lucid exposition again 
comes functional integral techniques which we sketch briefiy below. The four-point Green's 
function for the scattering of four massive scalar particles, interacting only gravitationally, 
is given by the time-ordered expression 

(O|r(0(xi)0(x2)0(x3)0(O)|O) = I [m]e^^^«G(xi, X3; g)Gix2, 0; g) + perms , (2.29) 

neglecting scalar loops, as they do not contribute to power behaved long range contributions. 
The propagators satisfy in the massless limit (it suffices to retain the massless limit as an 
approximation in the ultra-relativistic approximation). 
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- x'] g) = —6^^^ (x - x') , (2.30) 

where □ = g^^^'^d^{g^f'^g^'^dy) is the generally covariant Laplacian. If the metric is expanded 
about flat spacetime and the leading terms at high energies are retained, one obtains 

(-r^^'^d^d^ + V{x)) G{x, x') = 6'^^\x - x') . (2.31) 

In the leading approximation at high energy, the problem then reduces to potential theory. 
The potential is: 

V{x) = KE^h__{x) , (2.32) 

for the right mover and ( > +) for the left mover, where the light-cone coordinates are 

defined as usual by = z ± t (see in this context the papers by Kabat and Ortiz and 



Verlinde and Verlinde |^). One passes to the amputated proper one-particle propagator 
through the familiar scattering identities G = Gq + GqUGo = Gq + GqVG, where Go is the 
free one-particle propagator. Upon transformation to momentum space, one obtains 

(p'|S|j9) = (p'lV^l^p) = J rfxe-*f'-^/^r(x)$p(x) , (2.33) 

where "^pix) is the scattering solution which satisfies = \E'o + GqV'^ with appropriate 
scattering boundary conditions. Its WKB scattering solution is 

^p{x) = exp ^ |p ■ a; - ^ £^ F(x(r'))rfr'} , (2.34) 

where V{x) is given by Eq. (2.32), and the trajectories of the scattered constituents are 
given in light cone coordinates (+,tr, — ) by 

xi = {0,b,x~), X2 = (x+,0,0). (2.35) 

In Eq. (2.34), r is a coordinate that parameterizes the trajectory, r = x~^{x~) for right(left) 
movers, and 6 is a 2-vector representing the impact parameter. The on-shell scattering 
amplitude results from amputation of the Green's function, through use of Eqs. (2.29), 
(2.33), and (2.34), 
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A{E,b) = 2iE'^{Texp 



h2 



I K 



I 



(2.36) 



very much in parallel with Eq. (2.8). The energy momentum tensor in Eq. (2.36) now 
represents sources for the scattering of two relativistic particles in the so-called shock wave 
kinematics, 



The space-time coordinates in Eq. (2.29) have been simplified by the use of translational 
invariance and the form of the trajectories, Eq. (2.35) relevant to the scattering kinematics. 
The first term in Eq. (2.37) represents the right mover at impact parameter b and the second 
term represents the left mover at the origin in impact parameter space. 

A familiar form for the scattering amplitude is derived from Eq. (2.36) after quadratic 
functional integration to order G (k^) of the gravitational action. Upon Fourier transforma- 
tion to the transverse momentum representation, one obtains in the relativistic limit, the 
familiar eikonal formula with A = expiA, 



Equation (2.38) is the celebrated eikonal approximation corresponding to the functional 
average in Eq. (2.36); In lowest order the contribution to the eikonal is given by the two di- 
mensional graviton propagator i.e. A ^ {—2E'^G) InXb, which integrates to a form obtained 
in many references 



where a is the Regge trajectory for the Newtonian potential (a = 0(G|)) as s ^ oo), A 
is an infrared regulator for the 2-D effective massless propagator, and s and t are the usual 
Mandelstam variables. Eq. (2.39) is valid for large s and, near forward scattering, at small t. 
The real term in Eq. (2.39) is the purely classical scattering, while the pure phase (bracketed 
expression) is a contribution quantum mechanical in origin. The latter term corrects the 



T,, = E6'^'\x - a;i(r))5+e + E d'^'^x - x,{t))S;S; . 



(2.37) 




(2.38) 




(2.39) 
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ill-behaved growth of partial wave amplitudes at high s and their incumbent violation of 
unitarity 0. 

The general results Eqs. (2.36), (2.38), and (2.39)) have the desirable properties and 
features: 

• Unitarity in impact parameter space is retained at high values of s. Unitarity follows 
from reality of the exponential in Eq. (2.40) for elastic scattering and the positivity of 
its imaginary part for inelastic scattering for radiation of soft gravitons. See in this 
context the analysis of the CERN group in |]l9l . They establish a systematic expansion 



for the analogue of the static semi-classical results in G^s/lP'. 

The reader will see the close relation between the two contexts, ultra-relativistic and 
static. The ultra-relativistic theory is simply a 2-D image of the static 3-D effective 
theory. 

The reader will also recall that Eqs. (2.36), (2.38), and (2.39) together with the 
techniques used in their derivation were the basis of phenomenological models of high 
energy hadron scattering at least two decades ago [|T3|. Such techniques produced a 
geometric picture of high energy diffraction scattering of hadrons. 



III. RELATION TO STRING THEORY AND DISCUSSION 

The detailed calculations of the quantum corrections for example are presented only at 
the one-loop level, 0{h). It is obvious that in higher orders there will be difficulties. For 
example, divergent subgraphs by power counting will necessitate higher dimension polyno- 
mial valued counter terms coming from operators of dimension four and greater such as R^, 
where R is the Riemann tensor or any of its contractions. These divergent counterterms 
can then propagate into the skeletons in calculations analogous to those of Section III. It 
becomes impossible to proceed further in a theory that is divergent and non-renormalizable 
such as quantum gravity. Clearly some new definition of the theory is necessary. 
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The only well developed strategy is superstring theory in one of its various forms. We 
leave aside the vagaries of reconciling string theory to the phenomenology of particle physics 
and dwell on the question of a consistent theory of quantum gravity. Fortunately, there 
are one-loop calculations of on-shell scattering amplitudes in the literature p2| , p6| . There 
are even calculations beyond one loop that focus on unitarity which relate the higher loop 
graphs to tree graphs. The calculations are technical and involve the moduli space of higher 
genus Riemann surfaces, the general name given to certain complex manifolds that appear 
in string perturbation theory. The calculation of the one-loop amplitude by Montag and 
Weisberger is perhaps the most relevant, and we will quote and interpret their results 



for our purposes. 

In string theory another scale emerges; namely the Regge slope or string tension a' which 
plays two roles: It sets the scale of the coupling to and the mass scale for the higher mass 
states in the string spectrum of states. The field theory limit is then achieved by taking the 
a' —>■ limit. 

Higher-order graphs in string theory (of which there are fewer due to duality, e.g. one at 
the one loop level) involve integration over complex coordinates in contrast to the integration 
over momenta of the corresponding graphs in field theory. 

To order k^, the structure of the tree plus one-loop, on-shell, invariant scattering ampli- 
tude for massless string constituents i.e. dilatons or gravitons is 

M{s,t,u) = K[Ao{s,t,u) + cAi{a'; s,t,u)] , (3.1) 

where is a kinematical factor to ensure gauge invariance and proper mass-shell behavior, 
and c is a dimensionless constant. In Eq. (3.1), is the typical Veneziano-type amphtude, 

x_ 1 r(i - f s)r(i- ^t)r(i - ^^m) 
^' - r(i + f.)r(i + ft)r(i + f«) ' ^^'^^ 

while Ai{a'] s, t, u) is the one-loop string amplitude, the form of which is not essential for our 
purposes (it is given by an integral of certain world-sheet vertex operators over the moduli 
space of a torus). It suffices to note that in the low energy limit, the one- loop amplitude 
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can be expanded in powers of a' and the field theory hmit is recovered in the hmit a' 0, 
as long as \a's\, \a't\, \a'u\ <^ 1. The expansion is non-trivial due to subtleties in the region 
of integration in the moduli space which ensure a finite amplitude in the physical region 
of s, t, u. Nevertheless, when the dust settles, the result that emerges can be written as an 
expansion, where the first non-trivial term is 0{a'^) due to symmetric integration in moduli 
space is 



where A^'^'^ and A^"^"* are the first two terms in the expansion in a' . The first term A^^^ is 
the zero slope or field theory limit, and the next term is a local polynomial in the momenta 
at low energy, i.e. A^"^^ = 0(s^, t^, -u^). The field theory amplitude has a representation in 
terms of a Feynman parameter integral. 



where (f){s,t,l3) = /?i/?2'S — 13^13^1 is the usual polynomial in momentum invariants, d is the 
number of compactified dimensions (rf = 6), and there is to be a symmetrization with respect 
to the interchanges s — >■ t, s ^ m, etc. At low energy the massive states are completely 
decoupled and the field theory contributions completely dominate the amplitude. Hence it 
is seen that the field theory contributions overshadow the string corrections at low energies 
(large distances), as expected also from the effective action for the massless modes of the 
string (however, for a possible low-energy window on string-scale physics see p3[). 

Although the actual calculation presented here on the first non-trivial loop correction is 
clearly suggestive, we expect that higher order string calculations will support the general 
result of a gauge-invariant amplitude with polynomial corrections to the field theory result 
at low energies. At small distances, and s, t, and u large and comparable, the situation is 
quite different and string considerations are important and essential p^j35[] . 

We close this section with some further remarks concerning the relation of our work 
on the long range forces in quantum gravity with some related issues in gauge theories. 




(3.3) 




(3.4) 
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Namely, dynamical mechanisms of spontaneous symmetry breaking and mass generation 
render the Coulomb potential to be the only leading long range force in electrodynamics 
between charged objects. The finite fermion mass shields all long range forces at the quantum 
level in such theories. In non-abelian gauge theories such as QCD, confinement presumably 
renders all long range forces inoperable, once the hadron mass spectrum is achieved in a 
realistic calculational framework. 
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APPENDIX A: 

In this appendix we present some of the details of the calculations. The Feynman rules in 
de Donder gauge have been presented in the literature. We will present what is sufficient for 
the calculations at hand. The expansion of the metric Eqs. (2.5) and (2.10) give us all of the 
coupling to the matter sources needed. The one- and two-hoo{x) graviton vertices that couple 
to the matter sources are simply {—KiMj/2) and {iK^Mj/A), respectively. Furthermore, we 
need the triple graviton vertex. Actually, since we require the three-/ioo Green's function, 
we only need some special projections of the full triple vertex. One can show that the 
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momentum space transform of the configuration space three-point function involving only 
/loo's is given by the projection 

V3 — Hooooo ~ 2 (^000 M + ^0 ^ ij, 00 + fj, 0000) 

+ |(V^oo%% + V^%^oo + ^%oo^) (Al) 

where Vi^xnupaiP^Q, ~{p + q)) is the full triple graviton vertex in momentum space. The 
projected vertex is 

Vsip, q) = - J (p' + + (p + qf - 2(pogo - (Po + qof)) , (A2) 

These vertices are sufficient to compute the static potential to order G^. All contributions 
are relatively straightforward to compute except for the vacuum polarization of Fig. 3. All 
computations are performed in the harmonic gauge of Eq. (2.14). The contribution to the 
connected Green's function of Eqs. (2.17) and (2.18) from the graph in Fig. 2a is 

f2a = ^MiM2(Mi + M2) I dhdhdhGih, 0; ^2, 0; ^3, R) ■ (A3) 

where G(a;i, X2, ^3) is the three-point function of the temporal gravitons. Since all three 
gravitons are connected to static sources, the components of their momenta vanish, and 
we obtain that the three-point function in this case is the product of two propagators. 
The time integrals are trivial delta functions of po and go, times another delta function 
27r5(0) = T. Two of the three contributions to the final Fourier integral which transforms 
back to configuration space are divergent. However the divergence is a ultra-local divergence, 
which is irrelevant to the long-distance effects. The finite contribution is 

ha = ^M,M2(M, + M2)T / ' (A4) 

GM^M2 G(M, + M2) 
~ ~R R ■ ^^^^ 

Fig. 2b contributes 

f2b = \ (-y)' M,M2{M, + M2) J dtdhdhAit - h, R)A{t - R) . (A6) 
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A straightforward evaluation of the integral in Eq. (A. 6) yields the result 

Re J_l^ = + . (AT) 



This gives 



GM,M, GjM, + M,) 
/2. = ^T^ 2R • 



The sum of /2a and /2b gives rise to the post-Newtonian correction to the potential. 
Next, we consider the quantum corrections. First, the polarization contributes 

Pooa,n'^'''P,soo = (^) (- In(-g^)) . (A9) 

The vacuum polarization Fig. 3a receives contribution from both the graviton loop and the 
Fadeev-Popov-Feynman-Mandelstam ghost loop; this contribution has been computed more 



than two decades ago by Capper et al and Duff [0, and by 't Hooft and Veltman 
p8| . 't Hooft and Veltman obtain general results that can be transcribed to the long range 
potential that we are computing. The result can be written in terms of coordinate invariant 
quantities quadratic in the Ricci tensor and the curvature scalar. Donoghue has quoted 
their result in the form given below. We have not checked the complete calculation of the 
vacuum polarization from first principles. Recalling that the propagator is corrected as 

the vacuum polarization contribution to the quantum corrections follows. 
Fig. 3b is calculated easiest in configuration space, 

A' = i i^^jM.Ak J . (All) 

One notes that the square of the propagator is proportional to the derivative of a single 
propagator with respect to R. Recalling (2.19) and (2.20), we obtain 
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The quantum correction of the graph in Fig. 3c involves the same triple graviton vertex 
and calculates to 



/3c = -T 

Since 



32 



(27r) 



-ip-R 



1 



(A13) 



pO=0 



J 72^ 



{2nY p'^q^{p + qY 



p^ f d'^q 



pO=0 



12 J {2txY q^{p + qy 



pO=0 



12 



167r2 



ln(— p^) +p— independent^ , 



(A14) 
(A15) 



as can be shown in dimensional regularization, one obtains the advertised result, after dis- 
carding unimportant short-range divergences. It is useful to note that 



d^p 



hip — 



(A16) 



(27r)3 "'"^ 27ri?3 
Notice that other surviving zero mass non-gravitational contributions i.e. photons, mass- 
less neutrinos etc. will modify the results of this appendix, namely the vacuum polarization 
loop Eq. (A. 10), and produce additive contributions as expected. Massive fields will of 
course produce forces of finite range which do not survive at large distances. 

Notice also that the quantum corrections have a sign that increases the attractive grav- 
itational interaction between massive sources as the separation between constituents is de- 
creased. This sign has a physical origin. Namely the long range forces can be traced to real 
parts of amplitudes with positive definite absorptive parts coming from physical intermedi- 
ate states in the crossed channel which is the source of the logarithms; See Eq. (A. 15), for 
example. Unphysical polarizations in such intermediate states are cancelled by ghost loops 
as required by unitarity. Hence, positivity and the associated Schwarz inequality guarantee 
the definite sign of the three contributions of Fig. 3. While this is similar to quantum elec- 
trodynamics whose vacuum is dielectric in nature with dielectric constant greater than one; 
it is totally at variance with QCD. In QCD contributions purely real in origin reverse the 
sign of coupling constant renormalization and produce the inevitable asymptotic freedom 
and its associated paramagnetic vacuum. 
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